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Abstract 



A measurement of inclusive charged particle distributions in deep inelastic ep scatter- 
ing for 7*p centre-of-mass energies 75 < W < 175 GeV and momentum transfer squared 
10 < Q 2 < 160 GeV 2 from the ZEUS detector at HERA is presented. The differential 
charged particle rates in the 7*2? centre-of-mass system as a function of the scaled longi- 
tudinal momentum, xjr, and of the transverse momentum, pjf and <p\ 2 > , as a function 
of xf, W and Q 2 are given. Separate distributions are shown for events with (LRG) 
and without (NRG) a rapidity gap with respect to the proton direction. The data are 
compared with results from experiments at lower beam energies, with the naive quark 
parton model and with parton models including perturbative QCD corrections. The com- 
parison shows the importance of the higher order QCD processes. Significant differences 
of the inclusive charged particle rates between NRG and LRG events at the same W are 
observed. The value of <p^ 2 > for LRG events with a hadronic mass Mx, which excludes 
the forward produced baryonic system, is similar to the <Pt 2 > value observed in fixed 
target experiments at W fh Mx- 
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1 Introduction 



Inclusive particle distributions have been widely studied in deep inelastic scattering (DIS) 
and e + e~ annihilation to investigate the nature of the quark fragmentation and effects of higher 
order QCD processes. The formation of hadrons in DIS is a complicated process which cannot 
be fully calculated in the framework of perturbative QCD. In order to model this process it is 
convenient to distinguish two phases of the hadron formation. These correspond to a perturba- 
tive phase for QCD processes on the parton level followed by a non-perturbative fragmentation 
phase describing the confinement of the partons to observable hadrons. 

In this paper the charged hadron multiplicity distributions are analysed in the virtual-photon 
proton centre-of-mass system (7*^ cms), which corresponds to the centre-of-mass system of the 
produced hadronic final state with the invariant mass W . In the naive quark parton model 
(QPM) the virtual photon hits a quark in the proton and transfers a four momentum, q. The 
struck quark and the target remnant system each have an energy of W/2 in the j*p cms and 
move back-to-back with a 'velocity', which corresponds to a rapidity^] y max proportional to 
± In W. The outgoing quark and target remnant hadronise into multi-particle final states with 
limited p*, where p\ is the hadron momentum component transverse to the virtual photon 
direction as measured in the j*p cms. The width of the rapidity distributions of the produced 
hadrons is proportional to In W, while its height is approximately independent of W. From 
the measurement of jet profiles in DIS it is known that the width of a quark jet is typically 
two units of rapidity 0. At high values of W, the rapidity range populated by hadrons can 
be divided into three regions: the current jet region from (y max — 2) to y ma x, the region of the 
target remnant fragmentation from —y max to {—y max + 2) and a plateau region in between. 
When analysing hadron distributions as a function of the scaled longitudinal momentum in the 
7*p cms, xp, the current jet region defined above corresponds to the xp range xp > 0.05. If no 
QCD branching processes on the parton level are considered, the Xp and pi distributions for 
xp > 0.05 are predicted to scale in W. 

In fixed target DIS experiments J3], |], || effects of scale-breaking in the Xp distributions from 
QCD corrections, which are expected to soften the observed spectrum with increasing W, are 
small and could not be unambiguously identified. On the other hand, the mean square of pi, 
<pl 2 > , has been found to be very sensitive to higher order QCD effects |J. However, the 
details of the p* t spectra are also sensitive to non-perturbative fragmentation effects |BL 01, |S[]. 
With the high energies reached in ep collisions at HERA it is possible to extend the studies 
of xf and pi distributions to larger values of W, where the influence of perturbative QCD 
effects is expected to be much larger and the final state hadron distributions should reflect the 
dynamics of the subprocesses on the parton level. 

In a recent analysis the scaled momentum and charged multiplicity distributions of the hadronic 
final state were measured in the current region of the Breit frame as a function of the negative 
square of the four- momentum transfer, Q 2 , and the Bj or ken-scaling variable x ||. The charged 
particle spectra were observed to evolve with Q 2 in a way similar as in e + e~annihilation. In 
this paper we study inclusive charged hadron production as a function of xp and <p* t 2 > in the 
current region of the j*p cms frame. The objective of the analysis is to investigate the influence 
of perturbative QCD effects on the hadronic final state by studying the W dependence of these 
distributions in HERA ep collisions and in fixed target DIS data. The data are also compared 

1 The rapidity is measured with respect to the virtual photon in the 7*p cms. 



1 



with e + e~ results as well as with predictions of Monte Carlo programs. The comparison is 
also performed for a subclass of DIS events, which are characterised by a rapidity gap between 
the observed hadronic final state and the proton beam direction [TO] , and which are therefore 
candidates for diffractive scattering. 



2 The experiment 
2.1 HERA 

The data were collected during the 1993 running period using the ZEUS detector at the electron- 
proton collider HERA, where a 26.7 GeV electron beam and a 820 GeV proton beam were 
brought to collision providing an ep centre-of-mass energy of 296 GeV. 84 bunches were filled for 
each beam and in addition 10 electron and 6 proton bunches were left unpaired for background 
studies. An integrated luminosity of 0.55 pb _1 was collected. 



2.2 The ZEUS detector 



ZEUS is a multi-purpose magnetic detector which has been described elsewhere pTT| , [12|j . Here 
a brief description is given which concentrates on those parts of the detector relevant for the 
present analysis. 

Charged particles are tracked by the inner tracking detectors which operate in a magnetic field 
of 1.43 T provided by a thin superconducting solenoid. Immediately surrounding the beam pipe 



is the vertex detector (VXD) which consists of 120 radial cells, each with 12 sense wires [13 
The achieved resolution is 50 /im in the central region of a cell and 150 fim near the edges. 
Surrounding the VXD is the central tracking detector (CTD) which consists of 72 cylindrical 
drift chamber layers, organised into 9 'superlayers' [f[4}| . These superlayers alternate between 
those with wires parallel (axial) to the collision axis and those inclined at a small angle to give a 
stereo view. The hit efficiency of the CTD is greater than 95% and the resolution in transverse 
momentum for full length tracks is cr PT /p T = 0.005 px <S> 0.016 (pr in GeV), where ® means 
addition in quadrature. 

The solenoid is surrounded by a high resolution uranium-scintillator calorimeter (CAL), which 
is divided into three parts: forward^ (FCAL), barrel (BCAL) and rear (RCAL) [pL5|| . It covers 
99.7% of the solid angle. Holes of 20 x 20 cm 2 in the centre of FCAL and RCAL accommodate 
the HERA beam pipe. Each of the calorimeter parts is subdivided into towers which in turn 
are segmented longitudinally into electromagnetic (EMC) and hadronic (HAC) sections. These 
sections are further subdivided into cells, which are read out by two phototubes each. 

For measuring the luminosity as well as for tagging very small Q 2 processes, two lead-scintillator 
calorimeters are used ]12[ [16|]. Bremsstrahlung photons emerging from the electron-proton 
interaction point (IP) at angles 9 1 < 0.5 mrad with respect to the electron beam axis hit the 
photon calorimeter at 107 m from the IP. Electrons emitted from the IP at scattering angles 
less than 6 mrad and with energies between 20% and 90% of the nominal beam energy are 
deflected by beam magnets and hit the electron calorimeter placed 35 m from the IP. 

2 The proton beam direction is the forward +Z direction. 
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Two small lead-scintillator sandwich counters partially surround the beam-pipe at the rear of 
the RCAL. These counters were used to reject background produced by beam-gas interactions 
with the incoming proton beam and to measure the timing and longitudinal spread of both the 
proton and the electron beams of HERA. Two layers of scintillation counters mounted on either 
side of an iron veto wall, situated upstream of the detector, were also used to reject background 
particles. 



3 Data taking conditions 



The ZEUS trigger is organised in three levels [11| and reduces the input event rate from the 



bunch crossing rate of 10 MHz to 3-5 Hz. For DIS events, the first level trigger (FLT) requires 
at least one of three conditions for energy sums in the EMC calorimeter cells: the BCAL EMC 
energy exceeds 3.4 GeV; or the RCAL EMC energy (excluding the innermost towers surrounding 
the beam pipe) exceeds 2.0 GeV; or the RCAL EMC energy (including those towers) exceeds 
3.75 GeV. 

The second level trigger (SLT) rejects proton beam-gas events by using the event times measured 
in the rear calorimeter cells. The DIS trigger rate of the SLT is about one-tenth the FLT DIS 
trigger rate. The loss of DIS events at the SLT is negligible. 

The third level trigger (TLT) has the full event information available and applies physics-based 
filters. It requires tighter timing cuts to suppress beam-gas background further and also rejects 
beam halo muons and cosmic muons. The TLT selects DIS event candidates by calculating: 

5= Y, E i- (l-cos0i) > 20 GeV - 2£ 7 , 

i 

where Ei and 9i are the energy and the polar anglef] of the energy deposits in the calorimeter. 
The summation runs over all calorimeter cells. E 1 is the energy measured in the photon 
calorimeter of the luminosity monitor. For fully contained DIS events 5 ~ 2E e = 53.4 GeV, 
where E e is the energy of the incident electron. Photoproduction events have low values of 6 
compared to DIS events because the scattered electron escapes in the hole of the calorimeter 
which contains the beam pipe. 

For events with the scattered electron detected in the calorimeter, the trigger is essentially 
independent of the DIS hadronic final state. The trigger acceptance is greater than 97% for 
Q 2 > 10 GeV 2 and independent of Q 2 JXTJ. A total of about 7 ■ 10 6 events passed the TLT 
and was written to tape during the 1993 running period. 



4 Event kinematics 

In deep inelastic ep scattering events the incoming electron couples to a 7 or a Z (neutral 
current NC) or to a W + (charged current CC), which scatters off the proton. In the Q 2 range 
explored here, the contribution from W and Z exchange is negligible. The kinematic variables 
used to describe the inclusive DIS process are defined in Table |1| 

3 The proton beam direction is defined as the Z-axis in the HERA laboratory frame. 



3 



Variable 


Description 


1(1') 
P,M P 

u={P-q)/M p 
x = Q 2 /(2P-q) 
= Q 2 /(2M p u) 
y=(P-q)/(P-l) 
W 2 = {P + qf 

= Q 2 (l-x)/x + M 2 


Four-momentum of the incident (scattered) lepton 

Four-momentum of the proton and its mass 

Negative invariant mass squared of the exchanged virtual 

boson 

Energy of the exchanged boson in the proton rest frame 
Bjorken scaling variable 

Inelasticity parameter 

Invariant mass squared of the hadronic final state 



Table 1: Definition of the variables used to describe the kinematics of the inclusive DIS process 



The ZEUS detector is almost hermetic, allowing the kinematic variables Q 2 ,x and y to be 
reconstructed in a variety of ways using combinations of electron and hadronic system energies 
and angles. In the analysis presented here the double angle method (DA) was chosen, in which 



the scattered electron angle and the angle 7^ is used |]18fl . In the naive quark parton model 7^ 
corresponds to the angle of the scattered massless quark in the laboratory frame. The variable 
y is determined according to the Jacquet-Blondel method |L9| and is denoted by yjB- 



The four-momentum of the scattered electron needed to calculate the Lorentz boost to the 7*p 
cms frame, is reconstructed from its polar and azimuthal angle, 6 e ,<f) e . The scattered electron 
energy E' DA , used in the boost, is computed by the double angle method: 

E' DA = Q 2 DA /(2E e (l + cos6 e )) , (1) 

where E e is the energy of the incident electron and Q%a is given by: 

n 2 A F 2 sin 7H (l + cosfl e ) 

e ' sin lH + sin 9 e - sin ( 7h + 9 e ) ' 1 ' 



The variables xp and p* t describe the kinematics of the hadrons in the 7*7? cms: 

XF=PiM tmax \ = 2Pi/W , (3) 

where p* is the projection of the hadron momentum vector onto the direction of the virtual pho- 
ton and |p/* ma J is the maximum value of p\. The hadron momentum component perpendicular 
to the virtual photon axis is denoted by p\. 



5 Data selection 
5.1 Event selection 

The offline selection of DIS events was similar to that described in earlier publications 
(e.g. m Scattered electron candidates were selected by using the pattern of en- 

ergy deposition in the calorimeter. The electron identification algorithm was tuned for purity 
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rather than for efficiency. The purity is defined as the number of electrons generated and re- 
constructed in a bin divided by the total number of electron candidates measured in the bin. 
In studies with Monte Carlo DIS events and test beam data the purity was estimated to be 
> 96 % for E' DA > 10 GeV. 

The requirements for the final event selection were: 

• E' DA > 10 GeV, to minimise beam gas background contamination; 

• Qda > 10 GeV 2 ; 

• y e < 0.85, to reduce the photoproduction background, where y e is the scaling variable y 
as determined from the energy and polar angle of the scattered electron; 

• Vjb > 0.04, to guarantee sufficient accuracy for the DA reconstruction method; 

• 5 = — cos#j) > 35 GeV, where the sum runs over all calorimeter cells. For fully 
contained events 5 ~ 2E e = 53.4 GeV. This cut is used to remove photoproduction events 
and to control radiative corrections. 

Furthermore we required: 

• a primary vertex position, determined from VXD and CTD tracks, in the range 
-50 < Z vtx < 40 cm; 

• the impact point (X, Y) of the scattered electron in the RCAL to lie outside a square 
of 32 x 32 cm 2 centered on the beam axis, to ensure that the electron is fully contained 
within the detector and its position can be reconstructed with sufficient accuracy. 

After these cuts, the remaining photoproduction background was estimated to be ~ 1%. The 
contamination from beam-gas background was estimated to be below 0.5% as calculated from 
unpaired electron and proton bunches. Finally, QED Compton scattering events and residual 
cosmic and beam-related muons were rejected by algorithms, which identify this types of events 
by their pattern of energy deposits in the calorimeter cells. 

A total of 26100 events was selected by the above cuts. Of these events about 10% 0] contain 
a large rapidity gap in the hadronic final state. They are characterised by r] max < 1.5, where 
Vmax is the pseudorapidity of the most forward calorimeter cluster in the event, relative to the 
proton direction. The pseudorapidity is defined by r\ = — In (tan (0/2)) and a cluster is an 
isolated set of adjacent calorimeter cells with summed energy above 400 MeV. This sample 
is called the 'large-rapidity-gap' (LRG) event sample. The remaining events are denoted by 
'non-rapidity-gap' (NRG) events. The invariant mass of the hadronic final state excluding the 
scattered proton in the LRG events is calculated from the energy deposits measured in the 
calorimeter (excluding the electron cluster) by M x = \fY,had(E 2 ~Vx~Vy ~ Pz) ■ The values 
of px, Py and pz are the cell energies E projected on the axes of the HERA laboratory frame. 
The polar angles of these pseudovectors are calculated from the geometric centres of the cells 
and the primary event vertex position. The measured value of Mx is corrected to the hadron 
level as described in section |5|. 
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5.2 Track reconstruction and selection 



Tracks were recognised and fitted using two programs which were developed independently and 
follow different strategies for pattern recognition and track fitting. For the results shown in this 
paper the first approach is used and the second method was used for estimating the systematic 
error. 

In the first approach the track finding algorithm starts with hits in the outermost axial super- 
layers of the CTD. As the trajectory is followed inwards to the beam axis, more hits from the 
axial wires of the CDT and of the VXD are incorporated. A circle is fitted in the XY projection 
and is used for the pattern recognition in the stereo superlayer pattern. The momentum vector 
is determined in a 5-parameter helix fit. 



The second track finding program is based on the Kalman filtering technique |2"2fl. Seed tracks 



found in the outer layers of the CTD are extended inwards and points are added as wire layers 
of the CTD are crossed. The track parameters at each step are updated using the Kalman 
method. In the second step a Kalman fit to the points found in the pattern recognition phase 
is performed taking into account non-linear corrections to the measured drift time. Following 
the reconstructed CTD track inwards, CTD and VXD hits are associated with the track. The 
VXD track segments are merged with the CTD tracks using the Kalman filtering algorithm. 

Multiple Coulomb scattering in the beam pipe and in the walls of the VXD and CTD were 
taken into account in the evaluation of the covariance matrix. The vertex fit is performed with 
the fitted tracks using the perigee parameterisation j23fl . The vertex position is evaluated and 
the track parameters are calculated at the vertex. 

Only tracks which are associated with the primary vertex have been selected for this analysis. 
The tracks are required to have ptjab > 0.2 GeV and a polar angle in the HERA laboratory 
frame in the range of 25° < 9 < 155°. This is a region of the CTD, where the detector response 
and systematics are best understood. For tracks defined by these cuts the track reconstruction 
efficiency is ~ 95%. 

The scattered electron was removed from the track sample by rejecting those tracks which 
match the cluster in the calorimeter assigned to the scattered electron by the electron finding 
algorithm. Only tracks which reach at least the third superlayer and hence have a projected 
length in the plane perpendicular to the beam axis of more than 30 cm are kept to achieve 
the required transverse momentum resolution. For 9 > 150° the efficiency for identifying the 
scattered electron by matching the CTD tracks to energy deposits in the calorimeter decreases 
rapidly due to the limited acceptance and resolution of the CTD in the very rear part of the 
detector. Therefore the upper cut on 9 of the hadrons considered in the analysis was further 
tightened to 150°. 

Due to the cuts in 9 and pt,i a b the analysis in the 7*p cms is restricted to the range 10 < Q 2 < 
160 GeV 2 and 75 < W < 175 GeV, where the acceptance for charged hadrons is larger than 
60%. 

In Fig. [I] the distribution of the selected events in the Q 2 -x plane is shown. For comparison the 
kinematic region which has been investigated in fixed target experiments is also shown. 
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6 Acceptance correction 



6.1 Monte Carlo simulation 



The measured distributions are corrected for detector effects such as acceptance and resolution. 
For that purpose the hadronic final state from DIS was modelled using two different sets of 
Monte Carlo generators, the first for the description of the non-rapidity-gap events and the 
second to model the large-rapidity-gap events. 

Events from NRG DIS processes were generated using two alternative Monte Carlo models: a) 
the combination of the LEPTO 6.1 [|§ and the ARIADNE 4.0 Monte Carlo program |§] 



(CDMBGF) and b) LEPTO 6.1 with the option of combined matrix element and parton shower 



calculation (MEPS). The fragmentation was simulated using the LUND string model [27]] as 
implemented in JETSET |28| (see Table |J). 



Both models were interfaced to the program HERACLES p9| , which computes the electro- 
weak radiative corrections for DIS events. In the case of hard QED Bremsstrahlung the four- 
momentum vector of the virtual photon which probes the proton is significantly different from 
the virtual photon momentum reconstructed from the momenta of the incident and scattered 
lepton. In this case the xp and p\ distributions are also distorted and have to be corrected for 
this effect. In this analysis, however, the virtual photon momentum was reconstructed using 
the double angle method, which is insensitive to radiative effects. Events with hard QED initial 
state Bremsstrahlung photons [E hrems >, 7 GeV) are rejected by the cut on 5 > 35 GeV (see 
section 5.1). Monte Carlo calculations show that the QED radiative corrections are 5 — 10%. 

For both Monte Carlo simulations the MRSD'_ parameterisation of the parton densities in the 



proton was chosen [^0J, which gives a reasonable description of the structure function measured 
at HERA |31, 



The properties of LRG events are characteristic of diffractive interactions [Hj. Two Monte 



Carlo event samples have been used to model the hadronic final state of LRG events. The first 



was generated using the POMPYT Monte Carlo program ||33| , which is based on a factorisable 
model for high energy diffractive processes. Within the PYTHIA []34] framework, the incident 
proton emits a pomeron, whose constituents take part in a hard scattering process with the 
virtual photon or its constituents. The structure of the pomeron is assumed to be described 
by either a hard or a soft quark density function /(/?), where j3 denotes the fraction of the 
pomeron momentum carried by the quark. 

The second sample was generated following the model of Nikolaev and Zakharov (NZ) ||35|| , 



which was interfaced to the Lund fragmentation scheme In the NZ model it is assumed 
that the exchanged virtual photon fluctuates into a qq pair, which interacts with a colourless 
two-gluon system emitted by the incident proton. Both diffractive Monte Carlo samples were 
generated with default parameter settings. QED radiative processes were not simulated for 
these events. With the event selection cuts described in section 5, however, the QED radiative 
corrections are expected to be of the same size as for the NRG events. 

Event samples produced by the Monte Carlo generators marked in Table |2| by an asterisk were 
also processed by the ZEUS detector simulation program, which is based on GEANT 3.13 
P7| and which incorporates the detector and trigger simulation. Events fulfilling the trigger 
conditions were then passed through the standard ZEUS offline reconstruction program. 
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Acronym 



Description 



QPM 
CDM 
MEPS (*) 

CDMBGF (*) 

POMPYT (*) 



NZ 



Quark parton model + string fragmentation only 



Colour dipole model [25, 26 

Parton shower |M] matched to complete 0(a s ) matrix element 
calculation (ME) 

Colour dipole model combined with complete 0(a s ) matrix element 
calculation for the BGF process (ME) 

Model for diffractive DIS (assuming factorisation of the pomeron 
flux and the pomeron structure function) with 
a hard quark density function for the pomeron oc [/3(1 — j3)] 
or a soft quark density function for the pomeron oc [(1 — f3) 5 ] 
Model for diffractive DIS (non factorisable ansatz) l3l 



Table 2: Acronyms for the DIS models used in this report. For those generators marked 
by an asterisk, event samples have also been processed by the detector simulation and data 



reconstruction program. In all models the LUND string fragmentation model is used [27, 28 . 



The predicted r) max distribution for non-diffractive DIS events falls exponentially for i] max < 4, 
whereas for diffractive events this distribution is approximately flat. Calculations with the 
CDMBGF Monte Carlo model show that the fraction of non-diffractive DIS events with r] max < 
1.5 is about 5% fllQfl . The distributions for the LRG event sample defined by r] max < 1.5 
have been corrected with POMPYT and those for the NRG events have been corrected using 
the CDMBGF Monte Carlo program interfaced to HERACLES. Note that the results are not 
corrected for the selection inefficiency of the r\ max cut. 



6.2 Data correction procedure 



The measured hadron multiplicity distributions are distorted with respect to those of the true 
hadronic final state due to trigger biases, event and track selection cuts and the acceptance and 
resolution of the detector. The output of the trigger and detector simulation program together 
with the samples produced by the different event generators have been used to estimate the 
distortion of the distributions and to correct for them by multiplying the measured distributions 
by a correction function c(v ) in each bin of Q 2 and W, where v is the hadron variable under 
study and c{y ) is calculated as a bin-by-bin ratio: 



c(v) 



1 AN had (v) 



N, 



evt 



Av 



I 



1 AN had (v) 



gen 



N, 



evt 



Av 



(4) 



The subscripts gen and rec refer to the quantities as given by the event generator programs and 
the reconstructed quantities from the output of the detector simulation program, respectively. 
The number of events in a bin of Q 2 and W is denoted by N evt ; ANhad is the number of hadrons 
in a bin of v. The generated hadron distributions do not include the charged particle decay 
products of K°'s and A's and of weakly decaying particles with a lifetime > lCT 8 s. For the 
expression in the numerator events and hadrons are sorted in bins of the generated kinematic 
variables and for the denominator in bins of the reconstructed variables. In this way the 



S 



distributions have been corrected for losses of events and hadrons as well as for the effects of 
event migration, finite resolution and trigger biases. 



The bin size in the hadron variables v was chosen to be comparable with the estimated resolution 
in v and it was checked that the correction factor neither deviates by more than 40% from unity 



nor depends strongly on v [§8| . For models which adequately describe the data, the dependence 
of the correction factors on the model input was found to be small. The difference in c(v ) for 
different models was included in the systematic error. 



The mean square of pi (< p* t 2 >) was corrected by: 

h 

< Pt 2 >MC,rec 



< P f >=< k 2 > meas <Vt ^ >MC * en , (5) 



where <p* t 2 > me as is the mean value of pi 2 determined from the uncorrected data. The terms 
in the correction factor are defined as in equation |j. This method of correction is numerically 
more stable than the determination of < pi 2 > from the acceptance corrected pi 2 distributions. 

The following sources of systematic uncertainties were studied: 



• The model dependence of the correction factors c(v ) was estimated using two different 
models for the NRG and LRG event samples each. The CDMBGF and MEPS models were 
used to correct the NRG event sample and for the LRG event sample the POMPYT model 
with a hard quark density function (see Table 0) and the NZ model were used. The relative 
systematic error of 1/N evt ■ dNhad/dxF is ~ 3% and the one of l/N evt ■ dNhad/d <pl 2 > is 
-7%. 

• The analysis was done using two different strategies for track finding and vertex fitting 
as described in section 5.2. The difference of the corrected xf and pi distributions ob- 
tained with both programs is used as an estimate of the systematic error from the track 
reconstruction. The relative systematic error of 1/N evt ■ dN had /dxp is ~ 10% and the one 
of 1 /N evt ■ dN had /d <p* t 2 > is ~ 4%. 

• Systematic uncertainties in the determination of the four-momentum of the virtual photon 
may induce a systematic error in the hadron distributions measured as a function of xp 
and p* t . The size of this systematic error was estimated from Monte Carlo events by 
using the generated four-momentum of the virtual photon rather than the reconstructed 
four-momentum. The Lorentz transformation with the generated values was then used to 
calculate the momenta of the reconstructed final state particles in the 7*p cms and these 
values were compared to those obtained via the reconstructed virtual photon momentum. 
The relative systematic error of 1/N evt ■ dNhad/dxF is ~ 7% and the one of 1/N evt ■ 
dN had /d <p* t 2 > is ~ 5%. 

• The sensitivity of the measurements on the track selection criteria has been investigated. 
The cut in the polar angle of the tracks was varied between 20° and 33° and/or it was 
required that superlayer 5 instead of superlayer 3 has to be reached by the track. The 
requirement of a minimum hadron momentum transverse to the beam direction in the 
laboratory frame, pt,iah, was omitted. No significant changes in the results (< 1%) have 
been observed. 
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• The effect of a possible misestimation of the momentum resolution in the detector simu- 
lation program was studied by evaluating the correction function with a resolution of the 
measured transverse momentum artificially increased by 100%. The size of this effect on 
<pl 2 > and xf was smaller than 1%. 

The contributions of the above effects to the systematic error have been added in quadrature 
and are shown together with the statistical errors of the results in the tables and figures. 

The shape of the correction factors to be applied to the measured hadron distributions of xf 
and p* t as well as to < pi 2 > is shown in Fig. || separately for NRG and LRG events. The size 
of the correction for both event classes is very similar. 

7 Results 

7.1 xf and p* t distributions in NRG events 

First the xp and pi distributions of charged hadrons in NRG events are discussed. In Fig. [3]a 
the xp distribution at < W >= 120 GeV and < Q 2 >= 28 GeV 2 is compared with different 
models for hadron production in DIS. The Xp distribution falls steeply with increasing xp. The 
results from the HI experiment agree well with this measurements. The data agree with 
those models, in which higher order QCD processes are included, such as MEPS (solid line) 
and CDMBGF (dashed line), but not with the naive quark parton model (QPM) (dotted line). 

In Fig. |3|b the pi spectrum, which is integrated over xp > 0.05 for the study of the current 
jet fragmentation, is compared with the same model calculations. The QPM model predicts a 
much steeper pi distribution than the data show, whereas the MEPS model agrees well with 
the data. However, for closer investigation it is advantageous to take the mean square of pi, 

< p* t 2 >, a quantity which is more sensitive to the behaviour of the tail of the p* t 2 distribution. 

Figure ^c shows the <pl 2 > distribution as a function of xf for xf > 0.05. In any model, 
which allows for a transverse momentum of the partons, the rise of < pi 2 > with increasing 
xf is expected because a hadron with a higher value of xp carries also a larger fraction of the 
transverse momentum of the primary parton. Again the MEPS and CDMBGF models describe 
the data while the QPM strongly underestimates the value of <p* t 2 >. 

In Fig. ^a,c the results from the HI experiment are also shown |39fl . The differential hadron 
multiplicities measured by ZEUS are listed in Table ||]. 

7.2 xp and p* t spectra in LRG events 

The xf and p* t distributions from charged hadrons as well as <p* t 2 > as a function of xp are 
shown in Fig. (| separately for the samples of LRG and NRG events. The values for the LRG 
events are tabulated in Table The value of < W > is similar for both event samples, whereas 

< Q 2 > for the LRG events is lower by 30% than for the NRG events. The xp distribution for 
the LRG events is falling less steeply when compared to that of the NRG events. The LRG 
data in Fig. ^ are reasonably well described by the POMPYT (solid line) and the NZ (dashed 
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line) models for diffractive DIS with the t] max cut applied. The QPM prediction for the xp 
distribution of DIS events, shown by the dotted line in Fig. £|a, is slightly steeper than the xp 
distribution for LRG events. 

The p* t spectrum of LRG events is significantly less broad than that for the rest of the DIS 
events (Fig. |]b). This effect is highlighted in Fig. f|c. The mean values of pi 2 in events with a 
large rapidity gap are smaller than for the NRG events by a factor of 2-5. From a comparison 
with DIS model calculations with and without simulating QCD radiation processes, it is found 
that the <p* t 2 > values for LRG events resemble those for DIS events with only a small amount 
of gluon radiation. This observation is in good agreement with ZEUS results from the analysis 
of the energy flow ||21|| . However, < p* t 2 > in LRG events is somewhat larger than predicted by 
the QPM (see dotted line in Fig. ^c) , indicating that there is a non-zero contribution of higher 
order QCD processes in this class of events, too. This is confirmed by the observation of DIS 
events with a large rapidity gap which exhibit a two-jet structure p0[ . The model calculations 
for diffractive ep scattering slightly underestimate the measured values of < p\ 2 >. 

The inclusive distributions of LRG events have been found to have the properties of a diffractive 
interaction of a highly virtual photon with a proton ||10|| . Diffractive interactions in hadron- 
hadron reactions and photoproduction have been successfully described in the framework of 
Regge theory by the exchange of a pomeron [[41]. Several models have been developed to 



describe this reaction in terms of parton interactions (e.g. j43, [43]). In this context it is 
interesting to test the hypothesis that the diffractive DIS process can be viewed as the 'emission' 
of a pomeron from the proton, which carries the fraction x pom of the proton momentum, and a 
subsequent deep inelastic 7* pomeron scattering, which occurs at a higher value of x' = 



In this picture the relevant scale for the invariant mass of the hadronic final state should be 
given by Mj and not by W. In Fig. [5]a <p^ 2 > as a function of xp from the LRG events is 
compared with the results of a fixed target DIS experiment || , where the invariant mass of the 
total hadronic final state (<W>= 14GeV) is only slightly higher than the invariant mass of 
the hadronic final state observed in the LRG events (< Mx >= 8 GeV). The values of <p* 2 > 
and the xp distribution for both event samples agree reasonably well. This result supports the 
hypothesis that the transverse momentum space for the particle production is similar to DIS, 
where the scale of the invariant mass is given by Mx rather than by W. 



7.3 W and Q 2 dependence of xp and pi spectra 

In Fig. ||a the xp distribution from the NRG events is compared with that from e + e~ anni- 
hilation events on the Z° resonance [44], where the value of the e + e~ centre-of-mass energy 
is comparable to the value of W in the kinematic range analysed here. The differential rates 
for hadron production in e + e~ annihilation were divided by two so that they correspond to 
a single hemisphere and can be directly compared with the results from DIS. The differential 
hadron multiplicity distribution in DIS at HERA energies agrees with that observed in e + e~ 
collision events for xf^O.1. This confirms the approximate independence of the hadron forma- 
tion process from the type of the primary scattering objects, which most of the models assume 



27, 45 



The xp and <p\ 2 > distributions from this analysis are compared with those of DIS events at 
lower values of W [[K], [J7|. Since in fixed target experiments the DIS event sample has not been 
separated into NRG and LRG events, the NRG and LRG event samples have been combined 



11 



for the comparison. The xp and p* t distribution as well as <p* t 2 > as a function of xp for the 
NRG+LRG event sample are given in Table [5[ The distributions have been corrected using 
a combination of Monte Carlo event samples generated by the POMPYT and the CDMBGF 
Monte Carlo generator. The relative normalisation of the Monte Carlo samples has been fixed 
by fitting the sum of the reconstructed t] max distribution from the POMPYT and the CDMBGF 
Monte Carlo sample to the measured rj max distribution Iplfl . 

Figure |]b shows that the xp distribution becomes significantly softer with increasing W. The 
prediction of the QPM, where no scale breaking effects due to QCD radiation are included, 
almost agrees with the result from the fixed target experiments |7]. 46, 47 1 but is very different 
from the result at HERA energies (dotted line in Fig. The effects of scaling violation in 
the xp distributions of hadrons, which have been found to be small when measured in a limited 
interval of W and Q 2 0, |j, become evident when studied over a large range of W and Q 2 . 
Models in which higher order a s processes are considered (e.g. the MEPS model indicated by 
the full line in Fig. |6|b) agree reasonably with the ZEUS data. 

The mean value of p\ 2 as a function of xp is shown in Fig. [7| for < W >= 120 GeV (this analysis) 
and for < W >= 14 GeV from the EMC collaboration f46|. Comparing the results at low W 
and high W there is a strong increase of < p* t 2 > by a factor of about three over the whole 
range of xf > 0.05 going from W = 14 to 120 GeV. The comparison of the prediction from 
the QPM and the models including higher order QCD processes shows that QCD effects are 
much larger at HERA energies than at energies reached in fixed target experiments. 

For a further analysis of the W and Q 2 dependence, <pl 2 > was determined for two intervals 
in xp and four bins of W at an average value for Q 2 of 28 GeV 2 (Fig. ||a) and four bins of Q 2 
keeping W fixed at an average value of 120 GeV (Fig. ||b). The value of <pl 2 > increases both 
with W and with Q 2 . The results are tabulated in the Tables || and [7|. 

These results are compared with those from a fixed target experiment at lower energies |5|, 
[4"S|] . The rise of < p* t 2 > with W, which had been observed already in the fixed target DIS 
experiments, continues in the range of W seen at HERA. However, the Q 2 -dependence in these 
two ranges of W is different. There is a large overlap of the Q 2 intervals covered. At HERA 
energies a rise of <p* t 2 > with Q 2 is observed, while at low W almost no dependence on Q 2 was 
found H. 

The results from the ZEUS experiment and the fixed target experiment are compared with 
model calculations in Fig. |9|. The H^-dependence is reasonably described by the MEPS (solid 
line) and CDMBGF models (dashed line). Also the colour dipole model without including 
the BGF process (dotted line) qualitatively reproduces the W dependence of < p\ 2 > but 
overestimates the absolute value. The Q 2 -dependence is also described by the MEPS and 
CDMBGF model but not by the colour dipole model (CDM) alone. The colour dipole model 
simulates higher order gluon radiation processes but the BGF process is not considered. The 
Q 2 dependence of <p* t 2 > shows that it is necessary to include the explicit treatment of the 
BGF process as well in the simulation. 



8 Conclusions 

Measurements of differential charged hadron multiplicity distributions in DIS events have been 
presented in the centre-of-mass system of the virtual photon and the proton at a centre-of-mass 
energy of 296 GeV for 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV. 
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The transverse momentum, p* t , and xp distributions have been investigated separately for 
events with (LRG) and without a large rapidity gap (NRG) between the proton direction and 
the observed hadronic final state. In the whole range of xp > 0.05 the values of < p\ 2 > 
for NRG events are much larger than those for the LRG events. These results confirm that 
gluon radiation in LRG events is strongly suppressed as compared to 'standard' DIS events at 
comparable W. A comparison of the data with the prediction of the QPM shows, however, 
that some QCD radiation is present also in LRG events. 

The value of <p* t 2 > in the LRG events is similar to that observed in deep inelastic /ip scattering 
experiments on fixed targets at low W (< W>= 14 GeV). This indicates that the multi-particle 
production in LRG events is similar to that in DIS at a scale of the final state invariant mass 
W = M x , where M x is the invariant mass of the observed hadronic final state X, excluding 
the proton. 

The comparison of the Xp distributions in e + e~ annihilation and in DIS events confirms the 
hypothesis that the hadron formation process in the current jet region is approximately inde- 
pendent of the type of the primary interacting particles. 

The comparison of results presented here with those of DIS at low W from fixed target exper- 
iments allows a study of the development of QCD effects in the xf and pi distributions over a 
large range in W and Q 2 . A significant increase of <p* t 2 > with W is found. At HERA energies, 
the mean value of p* t 2 also rises with increasing Q 2 at fixed W. This can be understood in terms 
of the increase of the momentum space allowing the formation of more multi-jet events. 
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A 


Xp 


<Xp> 




1 dN had 




N ev t dx p 




0.03 




0.05 


0.04 




64.96 


± 0.78 


± 10.53 


0.05 




0.10 


0.07 




27.92 


± 0.34 


± 4.27 


0.10 




0.15 


0.12 




12.89 


± 0.23 


± 1.97 


0.15 




0.22 


0.18 




6.67 


± 0.13 


± 0.97 


0.22 




0.32 


0.27 




2.86 


± 0.07 


± 0.45 


0.32 




0.45 


0.38 




1.15 


± 0.04 


± 0.15 


0.45 




0.65 


0.52 




0.36 


± 0.02 


± 0.08 


0.65 




0.90 


0.73 




0.07 


± 0.006 


± 0.02 


B 


p\ GeV 


<p*> GeV 


Nevt dp; 


0.00 




0.10 


0.07 




1.05 


± 0.05 


± 0.18 


0.10 




0.20 


0.15 




2.84 


± 0.08 


± 0.51 


0.20 




0.40 


0.30 




3.89 


± 0.07 


± 0.60 


0.40 




0.60 


0.49 




3.34 


± 0.06 


± 0.47 


0.60 




0.80 


0.69 




2.24 


± 0.05 


± 0.33 


0.80 




1.20 


0.96 




1.10 


± 0.02 


± 0.18 


1.20 




1.50 


1.33 




0.47 


± 0.02 


± 0.07 


1.50 




2.00 


1.71 




0.23 


± 0.01 


± 0.03 


2.00 




2.75 


2.30 




0.08 


± 0.004 


± 0.01 


2.75 




3.50 


3.07 




0.03 


± 0.002 


± 0.004 


3.50 




5.00 


4.04 




0.01 


± 0.001 


± 0.001 


C 


Xp 


<Xp> 


<pf> GeV" 


0.05 




0.10 


0.07 




0.47 


± 0.01 


± 0.02 


0.10 




0.15 


0.12 




0.63 


± 0.02 


± 0.04 


0.15 




0.22 


0.18 




0.85 


± 0.04 


± 0.03 


0.22 




0.32 


0.27 




1.19 


± 0.06 


± 0.06 


0.32 




0.45 


0.38 




1.50 


± 0.09 


± 0.13 


0.45 




0.65 


0.52 




2.30 


± 0.19 


± 0.33 


0.65 




0.90 


0.73 




2.09 


± 0.33 


± 0.66 



Table 3: Differential multiplicities for charged hadrons as a function of A) xf and B) pi 
(xp > 0.05) and C) <p^ 2 > as a function of xp for DIS events with i] max > 1.5 (NRG) in the 
range of 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV. Statistical and systematic errors are 
given separately. 
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A 


Xp 


<x F > 


1 dN had 

Nevt dx F 


0.03 




0.05 


0.04 






26.97 


± 1.88 


± 2.55 


0.05 




0.10 


0.07 






19.70 


± 1.24 


± 0.94 


0.10 




0.15 


0.12 






12.16 


± 1.02 


± 1.51 


0.15 




0.22 


0.18 






6.49 


± 0.59 


± 1.01 


0.22 




0.32 


0.27 






4.02 


± 0.45 


± 0.82 


0.32 




0.45 


0.38 






1.41 


± 0.20 


± 0.46 


0.45 




0.65 


0.54 






0.68 


± 0.12 


± 0.10 


0.65 




0.90 


0.75 






0.31 


± 0.07 


± 0.07 


E 




p* t GeV 


<p*> GeV 


1 dN h d Q y-1 

Nevt dp* 


0.00 




0.10 


0.07 






1.36 


± 0.05 


± 0.50 


0.10 




0.20 


0.15 






4.31 


± 0.11 


± 1.11 


0.20 




0.40 


0.30 






4.91 


± 0.05 


± 0.41 


0.40 




0.60 


0.49 






3.44 


± 0.04 


± 0.34 


0.60 




0.80 


0.68 






1.45 


± 0.02 


± 0.55 


0.80 




1.20 


0.93 






0.64 


± 0.01 


± 0.045 


1.20 




2.00 


1.41 






0.09 


± 0.003 


± 0.02 


2.00 




5.00 


3.86 






0.001 


± 0.001 


± 0.004 


C 


Xf 


<Xf> 


<p* 2 > GeV" 


0.05 




0.10 


0.07 






0.19 


± 0.02 


± 0.01 


0.10 




0.15 


0.12 






0.24 


± 0.05 


± 0.03 


0.15 




0.22 


0.18 






0.35 


± 0.08 


± 0.04 


0.22 




0.32 


0.27 






0.41 


± 0.11 


± 0.16 


0.32 




0.45 


0.38 






0.38 


± 0.07 


± 0.11 


0.45 




0.65 


0.54 






0.50 


± 0.10 


± 0.12 


0.65 




0.90 


0.75 






0.37 


± 0.13 


± 0.60 



Table 4: Differential multiplicities for charged hadrons as a function of A) xp and B) p\ 
[xp > 0.05) and C) <p$ 2 > as a function of xp for DIS events with i] max < 1.5 (LRG) in the 
range of 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV. Statistical and systematic errors are 
given separately. 
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A 


Xp 


<Xp> 




1 dN had 




N ev t dx p 




0.03 




0.05 


0.04 




62.99 


± 0.75 


± 9.59 


0.05 




0.10 


0.07 




27.48 


± 0.33 


± 4.03 


0.10 




0.15 


0.12 




12.86 


± 0.22 


± 1.99 


0.15 




0.22 


0.18 




6.67 


± 0.13 


± 0.98 


0.22 




0.32 


0.27 




2.90 


± 0.07 


± 0.46 


0.32 




0.45 


0.38 




1.18 


± 0.04 


± 0.17 


0.45 




0.65 


0.52 




0.37 


± 0.02 


± 0.08 


0.65 




0.90 


0.73 




0.08 


± 0.007 


± 0.03 


B 


f t GeV 


<p*> GeV 


Nevt dp; 


0.00 




0.10 


0.07 




1.10 


± 0.05 


± 0.18 


0.10 




0.20 


0.15 




2.92 


± 0.08 


± 0.52 


0.20 




0.40 


0.30 




3.98 


± 0.07 


± 0.60 


0.40 




0.60 


0.49 




3.35 


± 0.06 


± 0.47 


0.60 




0.80 


0.69 




2.22 


± 0.05 


± 0.32 


0.80 




1.20 


0.96 




1.07 


± 0.02 


± 0.18 


1.20 




1.50 


1.33 




0.45 


± 0.02 


± 0.08 


1.50 




2.00 


1.71 




0.22 


± 0.01 


± 0.03 


2.00 




2.75 


2.30 




0.07 


± 0.004 


± 0.01 


2.75 




3.50 


3.07 




0.02 


± 0.002 


± 0.004 


3.50 




5.00 


4.04 




0.01 


± 0.001 


± 0.001 


C 


Xp 


<Xp> 


<p* 2 > GeV" 


0.05 




0.10 


0.07 




0.46 


± 0.01 


± 0.02 


0.10 




0.15 


0.12 




0.61 


± 0.02 


± 0.04 


0.15 




0.22 


0.18 




0.82 


± 0.04 


± 0.03 


0.22 




0.32 


0.27 




1.14 


± 0.06 


± 0.06 


0.32 




0.45 


0.38 




1.40 


± 0.08 


± 0.11 


0.45 




0.65 


0.52 




2.13 


± 0.18 


± 0.39 


0.65 




0.90 


0.73 




1.82 


± 0.28 


± 0.56 



Table 5: Differential multiplicities for charged hadrons as a function of A) xf and B) p\ 
(xp > 0.05) and C) < p* t 2 > as a function of xp for DIS events (combined NRG + LRG 
event sample) in the range of 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV. Statistical and 
systematic errors are given separately. 
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0.1 <x F < 0.2 


0.2 <x F < 0.4 


w 


GeV 


<W> GeV 


<p* t 2 > GeV 2 


<p* t 2 > GeV 2 


77 


- 95 


86 


0.59 ± 0.03 ± 0.07 


1.08 ± 0.08 ± 0.13 


95 


- 122 


108 


0.67 ± 0.03 ± 0.04 


1.19 ±0.08 ±0.10 


122 


- 141 


132 


0.74 ± 0.05 ± 0.04 


1.07 ± 0.09 ± 0.18 


141 


- 173 


157 


0.78 ± 0.05 ± 0.09 


1.44 ± 0.10 ± 0.12 



Table 6: <p* t 2 > as a function of W in two intervals of x F . Statistical and systematic errors 
are given separately. 





0.1 <x F < 0.2 


0.2 <x F < 0.4 


Q 2 


GeV" 


<Q 2 > GeV 2 


<p* t 2 > GeV 2 


<pf> GeV 2 


10 


- 20 


14 


0.59 ± 0.02 ± 0.03 


1.05 ± 0.05 ± 0.07 


20 


- 40 


28 


0.77 ± 0.05 ± 0.03 


1.27 ± 0.09 ± 0.04 


40 


- 80 


54 


0.86 ± 0.08 ± 0.10 


1.46 ± 0.16 ± 0.17 


80 


- 160 


110 


0.96 ± 0.11 ± 0.15 


2.12 ± 0.34 ± 0.77 



Table 7: <j»j 2 > as a function of Q 2 in two intervals of x F . Statistical and systematic errors 
are given separately. 
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Figure 1: Population of the Q 2 -x plane by the DIS events selected for this analysis. For the 
sake of clarity only 1/3 of the DIS event sample is shown in the scatter plot. Charged hadron 
distributions are investigated for 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV (dashed lines). 
The approximate kinematic region covered by the fixed target experiments is also indicated. 
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Figure 2: Inverse of the correction functions c(v) for a) the Xp and b) the p* t distribution in 
the range of 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV, which are used to correct the NRG 
event sample (full points) and the LRG event sample (triangles). The inverse of the correction 
function for < p* t 2 > as a function of xf in the same range of Q 2 and W is shown in c). 
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Figure 3: Differential charged hadron multiplicities for NRG DIS events normalised by the 
number of events as a function of a) xf and b) p* t for xp > 0.05. c) <p* t 2 > as a function of 
xf- For all plots the events are in a range 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV. The 
predictions of two DIS Monte Carlo models including QCD processes are shown: the MEPS 
model (solid curve) and the CDMBGF model (dashed curve). The prediction of the QPM 
is given by the dotted curve. The results of this analysis in a) and c) are also compared to 
measurements of the HI collaboration |39| . 
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Figure 4: Charged hardon distributions for 10 < Q 2 < 160 GeV 2 and 75 < W < 175 GeV 
(< Q 2 >= 28 GeV 2 and < W >= 120 GeV). a) The x F distribution, b) the p* distribution 
for xp > 0.05 and c) < p* t 2 > as a function of xp are presented separately for NRG and LRG 
events. In all three figures the curves represent the results of the following model predictions: 
solid curve: POMPYT with a hard pomeron structure function (see Table 0); dashed curve: 
model of Nikolaev and Zakharov; dotted curve: QPM. 
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Figure 5: Comparison of a) < p* t 2 > as a function of x F and b) the x F distribution for the 
LRG event sample (ZEUS) and DIS at low energy (EMC, < W >= 14GeV). The mean value 
of M x for the LRG event sample is < M x >= 8 GeV. 
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Figure 6: a) xf distribution from this analysis (NRG events) compared to results from e + e~ 
annihilation on the Z° resonance (W = 91 GeV) P3| . b) xp distribution from this analysis 
(NRG + LRG events) compared with results from fip DIS at < W >= 14 GeV |46| and at 
< W >= 18 GeV P7| . In Fig. 6b) the solid curve shows the prediction of the MEPS model 
calculation and the dotted curve that of the QPM at HERA energies. 
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Figure 7: <p* > as a function of xf from this analysis (NRG + LRG events) compared to 
results from \xp DIS at < W >= 14 GeV [TO. The curves show results from model calculations 
at HERA energy with the MEPS model (solid curve), the CDMBGF model (dashed curve) and 
the QPM (dotted curve). 
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Figure 8: <pjf 2 > in two intervals of function of a) W and b),c) Q 2 compared with 

results from \xp DIS experiments (EMC || and E665 The prediction of the MEPS Monte 

Carlo calculation is compared with the results of this analysis (solid curve) and of [H (dashed 
curve). 
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Figure 9: a) W and b) Q 2 dependence of <p* t 2 > from EMC || and ZEUS data compared with 
different model predictions for the hadron formation: MEPS (solid curve), CDMBGF (dashed 
curve) and CDM (dotted curve). For the explanation of acronyms see Table |2|. 
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